Fine Structure in the Decay of Deformed Proton Emitters: Non-adiabatic Approach 
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The coupled-channel Schrodinger equation with outgoing wave boundary conditions is employed 
to study the fine structure seen in the proton decay of deformed even-iV, oAd-Z rare earth nuclei 
131 Eu and 141 Ho. Experimental lifetimes and proton-decay branching ratios are reproduced. The 
comparison with the standard adiabatic theory is made. 
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Proton radioactivity has proven to be a very powerful 
tool to observe neutron-deficient nuclei and study their 
structure. Theoretically, proton radioactivity is an ex- 
cellent example of a simple three-dimensional quantum- 
mechanical tunneling problem. Indeed, in first order, 
it only involves a single proton moving through the 
Coulomb barrier of the daughter nucleus. In reality, the 
process of proton emission is more complicated since the 
perfect separation of the nuclear many-body wave func- 
tion into that of the proton and the daughter cannot 
be made, and - in addition - the decay is greatly influ- 
enced by nuclear structure effects such as configuration 
mixing. In spite of this, the first-order one-body picture 
works surprisingly well; it enables us to determine the 
angular momentum content of a resonance and the asso- 
ciated spectroscopic factor in many cases 0,^1 • Experi- 
mental and theoretical investigations of proton emitters 
are opening up a wealth of exciting physics associated 
with the coupling between bound states and extremely 
narrow resonances in the region of very low single-particle 
level density. One particular example of such a coupling, 
due to the Coriolis interaction, is discussed in this work. 

The last two years have seen an explosive number of 
exciting discoveries in this field, including new ground- 
state proton emitters and proton-decaying excited states 
|^-^|, and the first evidence of fine structure in proton 
decay [Q. The main focus of recent investigations has 
been on well-deformed systems which exhibit collective 
rotational motion; consequently, they are splendid lab- 
oratories for the interplay between proton emission and 
angular momentum. 

From a theoretical viewpoint, the understanding of 
proton emitters is a test of how well one can describe very 
narrow resonances. For spherical nuclei, there are many 
available theoretical methods, most of which give very 
similar and accurate results || . There have been several 
theoretical attempts to describe deformed proton emit- 
ters. These approaches can be divided into three groups. 
The first family of calculations is based on the 

reaction-theoretical framework of Kadmenskii and col- 
laborators poj . The second group is based on the theory 
of Gamow (resonant) states P,nfl-n3[. Finally, calcula- 



tions based on the time-dependent Schrodinger equation 
have recently become available All of these papers 
assume a strong coupling approximation. That is, the 
daughter nucleus is considered to be a perfect rotor with 
an infinitely large moment of inertia. Consequently, all 
the members of the ground-state rotational band are de- 
generate and the Coriolis coupling is ignored. Our work 
is the first attempt to go beyond these simplified assump- 
tions in the description of proton radioactivity. 

Our technique is based on the theory of Gamow states. 
More precisely, we solve the coupled-channel Schrodinger 
equation describing the motion of the proton in the de- 
formed average potential of a core (a daughter nucleus). 
It is assumed that the wave function of the proton is 
regular at the origin and asymptotes to a purely outgo- 
ing Coulomb wave. These boundary conditions result in 
complex-energy eigenstates For resonant states, the 
real part of the energy, Eq = Re(E), can be interpreted 
as the resonance's energy, while the imaginary part is 
proportional to the resonance's width, T — — 2Im(.E). 

Let us consider the Hamiltonian of the daughter-plus- 
proton system 

H = H d + H p + V, (1) 

where is the Hamiltonian of the daughter nucleus, H p 
is the proton Hamiltonian, and V represents the proton- 
daughter interaction. The total wave function, of the 
parent nucleus can be written in the weak-coupling form 
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In (|2|) u a (a=(Jdlpj P ) labels the channel quantum num- 
bers) is the cluster radial function representing the rela- 
tive radial motion of the proton and the daughter nucleus, 
yj p i p m p is the orbital-spin wave function of the proton, 
and &j d Md i s * ne wave function of the daughter nucleus. 
By definition, one has 
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d^J d M d - Ej d $j dMd - (3) 

In practice, the energies Ej d are taken from experiment 
or, if the data are not available, they are modeled theo- 
retically. Inserting (0) into the Schrodinger equation and 
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integrating over all coordinates except the radial variable 
r, one obtains the set of coupled equations for the cluster 
functions MM: 
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In Eq. (^) V a represents the average spherical potential 
of the proton in the state a, V£ a , is the off-diagonal cou- 
pling term, and Qj d is the energy of the relative motion 
of the proton and daughter nucleus in the state Jd- One 
obviously has Qj d = Qq — Ej d , where Qq is the Q p -value 
for the decay to the J]—Q + ground state. 

The method of coupled-channels described above has 
several advantages over the commonly used strong cou- 
pling formalism. First, excitations in the core may be 
included in a straightforward manner. This enables us to 
study the proton decay from the rotational bands of the 
parent nucleus to various rotational states of the daugh- 
ter. Furthermore, since the formalism is based on the 
laboratory-system description [Hamiltonian ([!]) is rota- 
tionally invariant and the wave function \1/ conserves an- 
gular momentum] , the Coriolis coupling is automatically 
included. 

The coupled equations (||) are solved in the complex 
energy plane. Asymptotically, the cluster wave function 
u Jdj i behaves like a purely outgoing Coulomb wave 

G i P i k J d r ) + iF i P ( k J d r ) with k J d = \Z 2 nQj d /h- In this 
work we assume that the average single-particle poten- 
tial is approximated by the sum of a Woods-Saxon (WS) 
potential, spin-orbit term, and the Coulomb potential. 
The axially deformed WS potential is defined accord- 
ing to Ref. JlTj . We employ the Chepurnov param- 
eterization Eq ]; it gives good agreement with proton 
single-particle energy levels as given in Ref. H]. The 
Chepurnov parameterization provides a reasonable com- 
promise between the Becchetti-Greenlees parameter set 
p9[ (excellent for the description of reaction aspects but 
slightly displacing the h%x/2, 99/2, and S1/2 proton shells) 
and the universal parameter set J2lj (excellent for the de- 
scription of structure properties of deformed rare earth 



nuclei |20| but having too large a radius to give a quan- 
titative description of the tunneling rate S). 

Since the resonance energy cannot be predicted with 
sufficient accuracy, following Refs. ||[|, the depth of the 
WS potential is adjusted to give the experimental Qq 
value. The deformed part of the spin-orbit interaction 
is neglected; we do not expect this to have a notice- 
able effect on the results [^2j . The off-diagonal coupling 
in (||) appears thanks to the non-spherical parts of WS 
and Coulomb potentials. The exact form of , can be 
found in Ref. |l6), Eq. (40) and Ref. §, Eq. (32). Here, 
the coupling potential is obtained by decomposing the 
WS potential into spherical multipolcs up to 12. 

We ensure that enough daughter states are considered 
for proper convergence. In practice, we must include 



some energetically forbidden states. These states do not 
directly contribute to the width, but do affect the solu- 
tion. Furthermore, we assume that the daughter nucleus 
is left in its ground-state rotational band and the defor- 
mation is unchanged during the decay process. To nor- 
malize the cluster radial functions, we use a method J23] 
which became known as "exterior complex scaling" |24| . 

The description of very narrow proton resonances is 
a challenging task due to dramatically different energy 
scales of Eq and T. Indeed, while the energies of single- 
proton resonances are of the order of 1 MeV, their widths 
can be as small as 10 -22 MeV. This calls for unprece- 
dented numerical accuracy. In this work, we apply the 
piece-wise perturbation method Q generalized to the 
coupled-channels case. The calculations are performed in 
extended precision arithmetic. The details of the numeri- 
cal procedure employed are given in Ref. p6| . As a check 
on the calculated widths, we also calculate the width from 
the probability current expression []l5| V = ^ Q T a , where 
the channel width is 



h 2 ul>* 



{r)u J a (r)-ui'(r)u J a *(r 



(5) 



The agreement between the two methods is always better 
than 0.1%. It should be noted that T is independent of r. 
For narrow and isolated resonances, one can approximate 
the exact E q. (pj ) by the R-matrix expression, as was done 
in Refs. $^\L?f 

One limit of Eq. (|I|) is the degenerate case in which 
Qj d = Qp f° r au values of Jd. This is the adiabatic ap- 
proximation discussed in Refs. |l6|,^7j. It is easy to check 
that in the adiabatic limit the set of new wave functions 
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with |if|<j p is also a solution of (Q) with an eigen- 
value Qq. For il=K=J, the wave function \&n = 

£ unn3p r lp( - r) y lpjpn represents the intrinsic single- 
particle Nilsson wave function with the angular momen- 
tum projection on the symmetry axis Q. As seen from 
Eq. (g), the strongly-coupled intrinsic state contains con- 
tributions from all the cluster wave functions correspond- 
ing to different core states. Another property of the adia- 
batic limit is the existence of solutions with J>f2. Since, 
as discussed by Tamura fl6|| , there is no dynamic coupling 
between the angular momentum of the proton and that of 
the daughter nucleus (the daughter nucleus is perfectly 
inert during the proton emission), there exist infinitely 
many solutions obtained by combining j p and Jd- Since 
the core states are degenerate, all the solutions with J>fl 
are degenerate as well. 

Let us discuss the results of our calculations. Since 
for the very proton-rich nuclei considered in this work 
practically nothing is known about their spectra, we pa- 
rameterize the ground-state band of the daughter nucleus 
as Ej = kJ(J + 1) and fix k to the experimental value 
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of E 2 + (or to the value taken from systematic trends). 
In the limit of infinite moment of inertia (re — > 0) one 
reaches the adiabatic limit. 

The presence of the rotational spectrum in the daugh- 
ter nucleus gives rise to rotational bands in the parent 
nucleus built upon the J—Q band-head. Figure [l] shows 
the calculated rotational band in 131 Eu built upon the 
J = 3/2 + level (associated with the [411]3/2 Nilsson or- 
bital). For the energy of the 2 + state in the daughter 
nucleus 130 Sm we took the experimental value 0). The 
J = 5/2 and 7/2 levels follow very closely the expected 
J( J + 1) spacing; the small deviations are due to the 
Coriolis coupling. 
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FIG. 1. Proton emission from the [411]3/2 rotational band 
in 131 Eu calculated in the non-adiabatic model. Half-lives for 
proton emission and branching ratios are indicated. 

The calculation of branching ratios (b.r.) is straight- 
forward in the non-adiabatic formalism. The partial 
width corresponding to the J—>Jd decay is given by 
T( J — > Jd) = J2j p i p ^J d j p i p ft * s secn tnat tne Proton- 
emission lifetimes and b.r. change with J. Of course, at 
the low energies shown in Fig. w, rotational states decay 
by emitting gamma radiation (r 7 3>r p ). (The competi- 
tion with proton emission is expected to take place at 
higher energies/spins.) The 3/2 + ^0 + decay is given by 
the very small d^/2 component. The large branching to 
the 2 + state is due to the dominating d 5 / 2 partial wave. 
(This also explains the fact that the 5/2 + level decays 
predominantly to the + ground state.) It is worth not- 
ing that the S1/2 component, not allowed in the adiabatic 
approach due to if-conservation, also contributes to the 
3/2+^2+ transition. Although K is not conserved in 
the non-adiabatic approach, one can decompose the wave 
function into states (^) with definite K . In most cases, 
we find that one i-T-component dominates the wave func- 
tion. Consequently, the Nilsson labeling convention can 
still be used. 

Transitions to excited daughter states, K^Jd, may 
also be approximated in the strong coupling framework 
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TABLE I. Half-lives and branching ratios (b.r.) to the 
jj~=2 + state for deformed proton resonances in 131 Eu and 
141 Ho calculated in the non-adiabatic (nad) and adiabatic 
(ad) formalism. The experimental values (shown in boldface) 
are taken from Refs. 17JJ5J. The energy of the 2 + state in 



130 Sm and 140 Dy was assumed to be 120 keV and 160 keV, 
respectively. 

Rpl. In this case the angular momentum conservation 
is guaranteed by the presence of the geometric factor 
(C'kk k-k) 2 ■ in addition, the Q p value is adjusted to 

As shown previously in Ref. |5j , at large deformations 
our calculations show a very small dependence on 02 and 
/?4. This is because the spherical decomposition of the 
corresponding Nilsson orbitals varies little in this regime, 
and there are no crossings between the levels of inter- 
est. The uncertainty due to the 02 value is usually much 
smaller than the experimental uncertainty in the proton 
energy. Table I shows predicted half-lives and b.r. for 
131 Eu and 141 Ho. For 131 Eu we take (3 2 = 0.32 and for 
141 Ho (P 2 = 0.29, /3 4 = -0.06) @. Spectroscopic factors 
have been estimated in the independent-quasiparticle pic- 
ture. Note that the l/(£! + 1/2) coefficient multiplying 
the BCS values of it 2 , assumed in Ref. is no longer 
present. 

Considering both the half- life and b.r., the ground 
state of 131 Eu is consistent with the [411]3/2 assignment. 
This result agrees with the analysis of Ref. and is at 
variance with the suggestion by Maglione et al. who as- 
signed the [413]5/2 orbital as a ground state of 131 Eu. 
The very small b.r. for the [413] 5/2 orbital results from 
the fact that both the 5/2+^0+ and 5/2+^2+ transi- 
tions go via the d 5 / 2 component which constitutes about 
4% of the wave function. On the other hand, for the yet- 
unobserved [532]5/2 state, the 5/2~^0 + transition goes 
via the tiny (0.1%) f$/2 wave, while the 5/2~^2 + decay 
is dominated by the / 7 / 2 component (17%) and the K- 
forbidden p 3 / 2 wave, which appears due to the Coriolis 
coupling. This results in the huge branching predicted 
for this state. 

For 141 Ho, based on calculated half- lives, we can assign 
the [523]7/2 level to the ground state and [411] 1/2 to the 
excited state. We note that these assignments are identi- 
cal to our previous assignment in Ref. M , although we are 
now using the non-adiabatic formalism and have changed 
the optical model parameters. These assignments also 
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agree with those proposed in Refs. ]3|jl3|] . 

For 140 Dy, the energy of the 2 + state is experimentally 
unknown, but it can be estimated from systematic trends. 
For instance, according to the N p N n scheme, one obtains 
the value i? 2 +=160keV pgj] , which was adopted in the 
calculations displayed in Table I. Figure || shows the ex- 
pected b.r. to the 2+ state as a function of E 2 + for both 
proton-emitting states in 141 Ho. For the [523]7/2 ground 
state, the predicted b.r. is still of the order of a few per- 
cent even at relatively large values of E 2 + , and this of- 
fers good prospects for its experimental observation. On 
the other hand, the b.r. for the isomeric [411] 1/2 state is 
lower by an order of magnitude. In 131 Eu the experimen- 
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FIG. 2. Predicted branching ratio for the proton transition 
from J^=7/2~ and 1/2+ states in 141 Ho to the 2+ state in 



Dy as a function of E 2 +- 
around 160 keV. 



The expected value of E 2 + is 



tal uncertainty in the Q p value leads to an uncertainty of 
+30% / —22% for all three levels. The branching ratios 
vary by less than 3%. 

In conclusion, we applied a non-adiabatic formalism, 
based on the coupled-channel Schrodinger equation with 
outgoing wave boundary conditions, to describe very nar- 
row proton resonances in deformed nuclei. The non- 
adiabatic model takes into account the fact that the 
daughter nucleus has a finite moment of inertia. Our 
calculations are consistent with the experimental data for 
the best deformed proton emitters known so far: 131 Eu 
and 141 Ho. As shown in Table I, the adiabatic approxi- 
mation gives rise to an underestimation of branching ra- 
tios sometimes by a factor 2-3. According to our predic- 
tions, there is a good chance to study experimentally the 
fine structure in the proton decay of 141 Ho. 
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